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In most magnetically-ordered iron pnictides, the magnetic moments lie in the FeAs planes, parallel 
to the modulation direction of the spin stripes. However, recent experiments in hole-doped iron 
pnictides have observed a reorientation of the magnetic moments from in-plane to out-of-plane. 
Interestingly, this reorientation is accompanied by a change in the magnetic ground state from a 
stripe antiferromagnet to a tetragonal non-uniform magnetic configuration. Motivated by these 
recent observations, here we investigate the origin of the spin anisotropy in iron pnictides using an 
itinerant microscopic electronic model that respects all the symmetry properties of a single FeAs 
plane. We find that the interplay between the spin-orbit coupling and the Hund’s rule coupling can 
account for the observed spin anisotropies, including the spin reorientation in hole-doped pnictides, 
without the need to invoke orbital or nematic order. Our calculations also reveal an asymmetry 
between the magnetic ground states of electron- and hole-doped compounds, with only the latter 
displaying tetragonal magnetic states. 


I. INTRODUCTION 

In the iron pnictides, unconventional superconductiv¬ 
ity appears in close proximity to a magnetic instability 
m- As a result, much of the research into these com¬ 
pounds has been devoted to understanding the magnetic 
properties of these systems [SHE]. Experimentally, the 
spin-density wave (SDW) magnetic order of most iron 
pnictides has orthorhombic ((7 2 ) symmetry and corre¬ 
sponds to stripes of parallel spins modulated either along 
the x direction (i.e. ordering vector Qi = (7T,0) and 
staggered magnetic order parameter Mi) or along the y 
direction (i.e. ordering vector Q 2 = (0, 7 r) and staggered 
magnetic order parameter M 2 ), in the coordinate system 
of the Fe square lattice 0EJ. Theoretically, this state has 
been described by a variety of approaches, from purely lo¬ 
calized Heisenberg spins [M3] to itinerant nesting-based 
scenarios [TOE121 to hybrid models mixing local moments 
and itinerant carriers [5] f23U27j . Common to nearly all 
these approaches is the assumption that the magnetic de¬ 
grees of freedom have an underlying 0(3) spin-rotational 
symmetry. From a phenomenological perspective, this 
implies that the magnetic free energy F mag depends only 
on the absolute value of the magnetic order parameters, 
i-e. Amag (Mf, M|) gS]. 

Despite the success of these approaches in describing 
many magnetic properties of the iron pnictides - such as 
the onset of a preemptive nematic transition |29] and the 
appearance of a tetragonal magnetic ground state [22] 
- there are important features that remain largely un¬ 
addressed. In particular, the 0(3) rotational symmetry 
of a free spin does not hold for a magnetic moment in a 
crystal. Instead, the symmetries of the underlying lattice 
induce anisotropies in spin space that may be significant 
[50] . Indeed, in most iron pnictides, the magnetic mo¬ 
ments are observed to point parallel to the modulation 


vector of the stripes, i.e. M, : || Qj [5J[B]. Attesting the 
significance of this spin anisotropy, a sizable spin gap of 
the order of 10 meV is also found at low temperatures 
deep in the magnetically ordered state [3TU34] . Inter¬ 
estingly, recent experiments in hole-doped iron pnictides 
have reported a spin reorientation near optimal doping, 
in which the direction of the magnetic moments flip from 
in-plane to out-of-plane [35H38] , Remarkably, this spin 
reorientation takes place in a region of the phase dia¬ 
gram in which the magnetic ground state changes from 
stripe/orthorhombic to tetragonal. 

Therefore, elucidating the origins of these spin 
anisotropies and their impact on the normal state prop¬ 
erties is paramount to advance our understanding of 
the iron pnictides. A natural candidate to account 
for these effects is the spin-orbit coupling (SOC) term 
AS • L which converts the lattice anisotropies 

into anisotropies in spin space. Recent ARPES measure¬ 
ments of the SOC A have reported values of the order of 
20 meV [43], which is not far from the typical magnetic 
energy scale of the problem (as extracted for instance 
from optical conductivity measurements 1441 [45] 1. To in¬ 
clude the SOC term in theoretical models, it is necessary 
to account for the puckering of the As atoms along the 
FeAs plane, which effectively doubles the unit cell of the 
Fe-only square lattice. In this paper, instead of working 
with the cumbersome ten-band model relevant for the 2- 
Fe unit cell, we consider a simpler low-energy microscopic 
model that respects all the symmetries of the FeAs plane 
and focuses only on the states near the Fermi level. Such 
a model, which relies on the smallness of the Fermi sur¬ 
face pockets of the iron pnictides, was previously derived 
by Cvetkovic and Vafek using rigorous group theoretical 
arguments [33] ■ Here, we show how the main ingredients 
of the model can be derived from a straightforward ex¬ 
pansion of the usual five-orbital model for the pnictides. 
By computing microscopically the magnetic free energy 
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in the paramagnetic state, we find the leading-order mag¬ 
netic anisotropic terms: 

SF = ai (M 2 X + Mly) 

+a 2 (M 2 y + M| x ) 

+ a 3 (M 2 Z + Mf,) . (1) 

The anisotropic coefficients are proportional not 
only to the square of the SOC term, A 2 , but also to the 
Hund’s rule coupling J. Evaluation of the coefficients 
reveals that a\ < 03 , 0:2 for most of the temperature¬ 
doping phase diagram, implying that the magnetic mo¬ 
ments have a general tendency to lie in the plane. In¬ 
terestingly, in the hole-doped side of the phase diagram, 
we find a small region in which 0:3 < 0 : 1 , 02 , indicating 
a spin reorientation from in-plane to out-of-plane. Both 
results are in qualitative agreement with the observations 
discussed above, providing evidence that the SOC term, 
with the aid of the Hund’s rule coupling, is sufficient to 
account for the magnetic anisotropies of the iron pnic- 
tides. This conclusion contrasts with previous proposals 
that orbital and/or nematic order are necessary to ex¬ 
plain the observed magnetic moment orientation [30] . 

For completeness, we also analyze the nature of the 
magnetic ground state across the phase diagram. We find 
a general tendency of electron-doped compounds to form 
an orthorhombic uniaxial (single-Q) stripe state (i.e. ei¬ 
ther (|Mi|) = 0 or (|M 2 |) = 0), whereas hole-doped com¬ 
pounds favor a tetragonal biaxial (double-Q) magnetic 
state (i.e. (|Mi|) = (|M 2 |)). Such an electron-hole asym¬ 
metry is also qualitatively consistent with experiments - 
and in particular with the recent observation that the 
spin reorientation takes place in a region of the phase di¬ 
agram in which the magnetic ground state is tetragonal. 

The paper is organized as follows: In section |TT] we in¬ 
troduce the low-energy microscopic model with the SOC 
term and the electronic interactions. Section ini is de¬ 
voted to the analysis of the coefficients of the free energy 
responsible for the magnetic anisotropy within leading- 
order. In section |IV| we refine the phase diagram by 
including fourth order contributions to the free energy 
that allow us to distinguish between stripe and tetrag¬ 
onal magnetic ground states. Concluding remarks are 
presented in section [V] Details of the calculations are 
included in four appendices. 


II. LOW-ENERGY MICROSCOPIC MODEL 


We start with a low-energy microscopic model that 
focuses only on the electronic states near the Fermi level, 
while respecting the symmetries of the FeAs plane. Such 
a model was originally derived in Ref. [39] using the 
symmetry properties of the non-symmorphic space group 
P4/nmm of a single FeAs plane. Here, we present an 
alternative derivation based on the typical 5-orbital tight- 


binding model used for the iron pnictides [461 : 

^0 = £ > lv (k) C jj,ka C i/,ka (^) 

k fiisa 

where k is the momentum, a is the spin, and y, v de¬ 
note one of the five Fe orbitals, xz, yz, x 2 — y 2 , xy, 
and 3 z 2 — r 2 . The matrix e(k) corresponds to the 
Fourier-transformed tight-binding dispersions involving 
up to fourth-nearest neighbor hoppings. Its explicit ex¬ 
pression is given in Appendix A. Note that this Hamil¬ 
tonian is based on the single-Fe square lattice (i.e. it 
refers to the “unfolded” Brillouin zone), and that the co¬ 
ordinate system is defined such that k x and k y are par¬ 
allel to the nearest-neighbor Fe atoms directions. The 
actual crystallographic unit cell contains two Fe atoms 
due to the puckering of the As atoms, resulting in the so- 
called “folded” Brillouin zone, described by the coordi¬ 
nates K x ,K y (see Fig. [I]). Note that the two coordinate 
systems are related by: 


A ; — k x T ky 
Ky — k x T . 


(3) 


where the momentum in the unfolded zone is measured 
in units of its inverse lattice constant 1 /a, whereas the 
momentum in the folded zone is measured in units of the 
its inverse lattice constant 1 / (i/2a). 

The key properties that allow us to derive a simpler 
low-energy model are the facts that the Fermi surface 
pockets are small and that the orbitals that mostly con¬ 
tribute to the Fermi surface are xz, yz, xy. In particular, 
the idea is to start at the high-symmetry points of the 
unfolded Brillouin zone (namely, T = (0,0), X = ( 7 r, 0 ), 
and Y = ( 0 , 7 r)), where the band states are pure orbital 
states, and perform an expansion of the corresponding 
matrix elements £^ (k) for small momentum. Note that, 
to focus on a general and analytically tractable model, 
we follow Ref. [35] and ignore the states near the (n, ir) 
point of the unfolded Brillouin zone. While it is true that 
some iron pnictides display a hole-pocket with sy-orbital 
character centered at this point, this pocket is not usu¬ 
ally present for all values of k z , and is absent in many 
of the iron-based materials with a single FeAs plane per 
unit cell. Correspondingly we consider in this work the 
doping range in which this pocket lies below the Fermi 
level. 

Consider first the F point; the two states closest to 
the Fermi level are the xz and yz orbitals, which form 
a degenerate doublet in the absence of SOC. Thus, for 
small k, we define the spinor: 


Vr.k 


/ CyzM \ 

c xz, kt I 
c yzMl I 

V Cxz,\<l\, / 


(4) 


Projecting e^ (k) on this sub-space and expanding for 
small k then yields the 4x4 Hamiltonian: 


#o,r = V’r.k^r (k) V> r ,k (5) 

k 
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with 


M k ) = 

(^+ 2 ^+ b ( k l~ k l) ick^ky 

\ 4 ck^ky e r +2 ^.-b(kl-k 2 y ) 


( 6 ) 


where er° is a Pauli matrix acting on spin space. The co¬ 
efficients er, mr> b, and c can be obtained directly from 
the tight-binding parameters (see Appendix A). Note, 
however, that as we move away from the high-symmetry 
points of the Brillouin zone, other orbitals start to con¬ 
tribute to the electronic states. Consequently, the coeffi¬ 
cients of the expansion (as derived in Appendix [AJ will 
be slightly renormalized by the hybridization with the or¬ 
bitals not included in the expansion, although the form 
of the expansion remains invariant. To account for this 
issue, we can consider the coefficients to be free param¬ 
eters that can be fit directly to the first-principle band 
dispersions. 

Near the X point, the low-energy states correspond to 
the orbitals yz and xy. Defining the spinor: 


Vbc.k+Qi 


Cyz, k+Qit 
Ccy.k+Qif 
Cj/z,k+QiJ. 
Cxy, k+Qit 


(7) 


and expanding the projected e (k) near Qi = (w, 0) 
yields: 


Ho,X — ^A'.k+Qi b'X (k + Qi) V’.Y.k+Qi (®) 

k 


with 


h. v(k + Qi) = 

/ e l+ 2 ^+ a i( k l~ k l) -*«x(k) 

\ *«x(k) £ 3+2 2^TJ+a 3 (fcJ-fcy) 


and: 


(9) 


u.Y(k) = 2vky + 2p 1 k y (ky + 3kl) - 2p 2 k v (kl - k y ) . (10) 

Similarly, near the Y point, the low-energy states in¬ 
volve the orbitals xz and xy: 


IpY, k+Q 2 


f Ca, Zj k+Q2t 
Qcy,k+Q 2 t 
C ai 2 :,k+Q 2 i 
Cxy, k+Q 2 t 


(ii) 


Projecting and expanding e^ (k) near Q 2 = (0, it) gives: 



FIG. 1. (Color online) (Left) Illustration of the 2-Fe (dotted 
line) and 1-Fe (solid line) unit cells. The black dots denote 
iron atoms, while the pnictogens form two sublattices, one 
above the iron-plane (dark blue) and one below the iron-plane 
(light blue). (Right) Brillouin zones corresponding to the 1-Fe 
and 2-Fe unit cells. The dotted line is the “folded” Brillouin 
zone, corresponding to the 2-Fe unit cell, while the solid line 
is the “unfolded” Brillouin zone, corresponding to the 1-Fe 
unit cell. 


and: 

uy(k) = -2vk x -2p 1 k x (kl+3kl)-2p 2 k x (kl-kl). (14) 

Having established the low-energy states in the un¬ 
folded Brillouin zone (i.e. the one referring to the 1-Fe 
unit cell), it is now straightforward to fold the states 
into the 2-Fe unit cell (see Fig. [l]). Despite working in 
the folded Brillouin zone, described by the coordinates 
K x ,K y , we will still make use of the coordinates k x ,k y 
of the unfolded zone. From Eq. (pi), we find that upon 
folding, both momenta Qi = (71-, 0) and Q 2 = (0, 7 r) are 
identified with the same momentum Q m = ( 7 r, 7 r). It 
is straightforward to show that the spinors X and Y 
now combine to form two new degenerate doublets at 
the M = ( 7 r, 7 r) point of the folded zone: 


C xy,k+Q 2 t 
C xy,k+Qi| 
c xy, k+Q 2 | 
C xy, k+QiJ. 

(15) 

Hereafter, we will consider the momentum of any 
spinor as measured relative to the high-symmetry points, 
as appropriate. Then, the non-interacting Hamiltonian 
becomes: 

^o = E^k[ ff o(k)-Mll]^ k , (16) 

k 

where we defined the enlarged spinor: 


V>M 1 ,k+Q M = 


/ Ca;z,k+Qst 
c yz, k+Qif 
Cxz,k+Q 2 4. 
V c yz ,k+Qi 1 


;^M 3 ,k+Q„ = 


Ho,Y — y] V , y ! k+Q 2 ^Y (k + Q 2 ) Vv,k+Q 2 P^) 

k 


with 


hy( k + Q 2 ) — 

/ £ i+ 2 2^r^ a i { k l~ k D -wv(k) 

\ ivy (k) e 3+ 2 !^-“3(^-^) 


u 0 (13) 


/ 'tPYM \ 
ikk = V'A.k 
V ^r,k ) 

and the Hamiltonian matrix: 

/Mk) 0 0 \ 

TLo(k) = 0 h x ( k) 0 

V 0 0 h T (k)J 


(17) 


(18) 
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where y is the chemical potential and 1 is the identity 
matrix. We note that this model has the same properties 
of the Hamiltonian derived by Cvetkovic and Vafek in 
Ref. [391 combining a k • p expansion and the symmetry 
properties of the P4/nmm space group (note, however, 
that the definition of the spinors X and Y are switched 
in Ref. [32] with respect to the notation adopted here). 
In the group-theory language, the spinor ipr belongs to 
the two-dimensional E g representation of P4/nmm near 
the r point, whereas ^Mi and ipM 3 belong to the two- 
dimensional Em 1 and Em 3 representations of PA/nmm 
near the M point. Hereafter, we will use for the coeffi¬ 
cients of the Hamiltonian the parameters given by Table 
IX in Ref. [39]. Those were obtained by direct fitting of 
the band dispersions to first-principle calculations. The 
resulting band dispersions, as well as the Fermi surface, 
are shown in Fig. [2] Note also that this low-energy 
model is fundamentally different than two-orbital mod¬ 
els that restrict the Hamiltonian to the subspace of the xz 
and yz orbitals. Our model, derived from the five-orbital 
tight-binding model as shown in Appendix A, not only 
obeys all the symmetries imposed by the P4/nmm space 
group, but it also contains information about all the or¬ 
bitals that contribute to the Fermi surface, including the 
xy orbital. 

Besides the band dispersions, the non-interacting 
Hamiltonian must also contain the SOC term AS • L, 
with S denoting the spin angular momentum operator 
and L, the orbital angular momentum operator. Note 
that this atomic-like term preserves the Kramers degen¬ 
eracy of each state. To proceed, we project this term 
from the L = 2 cubic harmonic basis to the orbital basis 
(see Appendix B for more details). At the T point, we 
obtain an admixture of the xz and yz orbitals: 

\ H { iC lz.W><fJ C xzMJ + hx -) = 

k a/3 

^X^r.kO^ ® O^r.k ( 19 ) 

k 

where, in the last step, we used the definition of the 
spinors. At the M point, we obtain the admixture of 
xz/yz and xy orbitals: 




FIG. 2. (upper panel) Cut of the low-energy band dispersion 
from the T = (0, 0) to the M = (n, tt) point of the folded 
Brillouin zone with parameters fit to the tight-binding model 
of Ref. [2?]. The corresponding Fermi surface is shown in the 
lower panel. 

with r ± = ,] ( t x ± it v ). Therefore, the SOC becomes: 

'Hsoc = ^k-Hsoc(k)4' k i (22) 


£(‘ 


'xz,\aa. u 


kcn/3 


E [*4-,: 


k+Q 2 


■ ® <j x ) ip 


A'.k+Q, 


+ h.c. 


( 20 ) 


as well as: 


A 

2 


( ic ly,k* a lp c yzW + hx -) 


ko:/3 


J2 [*4,k+Q 2 ( r ® ° v ) Vbck+Q! + h-c.] (21) 


with: 

( 0 h s °c (k) o \ 

ffsoc(k) = (/4 oc ( k)) f 0 0 (23) 

V 0 0 /if° c (k)/ 

such that: 

/ip OC (k) = ?, A (r v ® a z ) , (24) 

h^f c ( k) = ^A (r + (g) <j x + t~ <g> cr y ) , (25) 

in agreement with the group-theoretical arguments of 

Ref. [32]. 
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The interacting part of this low-energy model is rather 
complex, involving 30 different possible biquadratic 
terms in the fermionic operators. Here, we will focus on 
the interactions coupling the T and the M points, since 
those are the ones that will be relevant for the calcula¬ 
tion of the magnetic action in the next section. Defining 
t 1 " 3 = ^(t° ± r z ), the interacting terms coupling the T 
and M points are written as (see Ref. [32]): 


which, transformed to the spinor representation, yields: 


-Hsdw = Mi • [V’r.k ( f 1 ® °0 ^xm+Qi + hx ' 

k 

+M 2 • [l4,k ® or) ^y,k+q 2 + h c - 

k 


(28) 


Note that, in the language of Ref. 59], the SDW 
Hamiltonian transforms under the Em 4 two-dimensional 
irreducible representation of PA/nmm. 


Pint 




V13 (v4<r( k ) T k) + h.C.) 


+V 13 (^Ycr (k)f 3 ?/>rv (k) + h.c.) 
+*’i 5 (^ CT (k)f 3 ^r CT (k) + h.c.) 
+V 15 (V , y CT ( k )' r_ ^r I r(k) + h.c.^ 
+«i 7 (V , x CT ( k )'r + ^r CT (k) + h.c.) 
+^i 7 (V’ivO^Vr^k) + h.c.) 
+ui 9 (V'L( k ) f V^( k ) + h.c.) 
+U19 (V’y<r( k ) T+ V>rv(k) +h.c.) , 


2 


(26) 


Note that all terms are diagonal in spin space. In terms 
of the more usual multi-orbital Hubbard model with on¬ 
site interactions, the first three coefficients originate from 
the Hund’s rule coupling, U 13 = V\ 3 = Vn = J, while 
the last one arises from the intra-orbital Hubbard term, 
U 19 = U /2 [39] . Here, we are not interested in which in¬ 
teractions will drive the SDW transition. Rather, we will 
assume a nearby SDW instability and compute how the 
interplay between these interactions and the SOC affect 
the magnetic action. 

Finally, in order to be able to derive the magnetic ac¬ 
tion in the next sections, we need also to establish how 
the magnetic order parameters Mi and M 2 , correspond¬ 
ing to ( 7 r, 0) and (0,7 r) order in the unfolded zone, cou¬ 
ple to the low-energy electronic states. Here, we will 
consider only intra-orbital magnetism. Indeed, previous 
Hartree-Fock investigations of the five-orbital Hubbard 
model have shown that the dominant contributions to 
the magnetic instability arise from intra-orbital couplings 
[581 [591 . Therefore, the SDW vertices become: 

-Hsdw = Mi ■ ( c lz,ka (To ‘P c yz ,k+QiP + hx -) 

k ctfi 

+M 2 ■ ( c L,ka <T «^c x2;k+Q2 ^ + h.c.) (27) 

k af3 


III. ANISOTROPIC MAGNETIC FREE 
ENERGY 

To understand the origin of the magnetic anisotropies, 
we first review the group theoretical arguments of Ref. 
[39] , In the absence of SOC, all the components of the 
magnetic order parameters belong to the same irreducible 
representation , as shown above, and the free energy 
depends only on the invariant form M| + M|. How¬ 
ever, with the introduction of SOC, the spin and orbital 
degrees of freedom are no longer independent. Conse¬ 
quently, the components of M ; must belong to differ¬ 
ent irreducible representations of P4/nmm, if one en¬ 
forces the combination of the spin and orbital parts of 
the magnetic order parameter to still transform under 
E Mi . As a result, the individual components of Mi and 
M 2 transform according to the following two-dimensional 
irreducible representations |39] 



Because these components belong to different irre¬ 
ducible representations, they will, in general, have differ¬ 
ent transition temperatures. Therefore, the free energy 
must acquire the leading-order anisotropic terms: 


SF = a.\ (. Ml x + MfJ 

+a 2 (Mi \ y + 

+a 3 (. M*, + Ml,) . (30) 

The smallest cci coefficient determines which type of mag¬ 
netic order condenses first. In Fig. [3] we show separately 
the real-space spin configurations corresponding to the 
components of Mi and M 2 associated with each coef¬ 
ficient a*. Specifically, if a 3 < a 2 ,a 3 , then M, points 
parallel to the ordering vector Q, (of the unfolded zone); 
if a 2 < on, 0 : 3 , M, still points in-plane, but perpendicu¬ 
lar to the ordering vector Q;. Finally, if a 3 < a 3 ,a 2 , M, 
points out-of-plane. Note that this analysis does not re¬ 
veal whether only either Mi or M 2 condense, or if both 
condense simultaneously. To establish the actual ground 
state, it is necessary to go to higher order in the free 
energy. We will come back to this point in Section |IV| 
Note that the spin-anisotropic terms preserve the tetrag¬ 
onal symmetry of the system. 
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FIG. 3. Sketch of the different uniaxial (i.e. stripe-like) mag¬ 
netic configurations corresponding to the different anisotropic 
terms in the magnetic free energy (301 with coefficients a\, 
« 2 , and 0 : 3 . 



FIG. 4. Schematic representation of the two distinct type of 
Feynman diagrams at 0(v, A 2 ): two-loop diagrams (a) and 
one-loop diagrams (b). Only the one-loop diagrams con¬ 
tribute to the anisotropic terms. 


Here, our goal is to evaluate microscopically the at co¬ 
efficients using the model of the previous section. Within 
the non-interacting part of the model, we find that even 
the presence of spin-orbit coupling does not introduce 
magnetic anisotropies. The reason is that the model ef¬ 
fectively has an enlarged P4/mnm®P4/n?nm symmetry, 
since the states at T and the states at M are treated in¬ 
dependently. Of course, the fact that these states are 
connected in realistic tight-binding models ensures that 
some level of spin-anisotropy will be introduced at the 
non-interacting level. Such an effect will be likely a high- 
energy effect, as it involves states away from the Fermi 


level m- Here, instead, we focus on the low-energy con¬ 
tributions to the spin anisotropy. Consequently, they 
must come from interactions - particularly, from the in¬ 
teraction terms that couple the states at F and at M, 
and therefore remove the enlarged PA/nmm PA/nmm 
symmetry. These are precisely the terms listed in Eq. 


We proceed with a straightforward diagrammatic ap¬ 
proach by dressing the non-interacting particle-hole bub¬ 
ble with the SOC term A in Eq. (231 and with the in¬ 
teractions Vi in Eq. (26). The SDW vertices coupling 
the magnetic order parameters to the non-interacting 
Green’s functions are those derived in Eq. (28). Both 
A and v 2 are treated perturbatively to leading order. Be¬ 
cause terms of the order 0(A) are forbidden by symmetry, 
we consider the diagrams of the orders 0 ( A 2 ) and G(vi). 

To order G{vt ), there are two distinct types of 
interaction-dressed diagrams, as depicted in Fig. [4j On 
top of that, to order 0 { A 2 ), each of the two diagrams can 
be dressed by a pair of SOC legs in eight different ways. 
Because symmetry requirements forbid terms that cou¬ 
ple directly M 4 and M 2 at the quadratic level, the pair 
of SOC legs must correspond to the same SOC term, i.e. 
either /ip OC (k) or /i|p c ( k) in Eq. (23). Explicit calcula¬ 
tion of the traces over the Pauli-matrices reveals also that 
the only combinations of SOC legs that yield anisotropic 
magnetic terms are those in which one SOC leg appears 
in the upper-right (lower-right) part of the diagram and 
the other SOC leg appears in the lower-left (upper-left) 
part of the diagram. 

We find that all two-loop diagrams (i.e. those rep¬ 
resented in Fig. So vanish, and therefore do not con¬ 
tribute to the magnetic anisotropy term (30). We show 
this explicitly in Appendix 0 Therefore, all that is left 
is to compute the one-loop diagram represented in Fig. 
[4Jd. The calculation is tedious but straightforward. To 
illustrate it, consider the interaction U 13 . The one-loop 
diagrams contributing to the anisotropic terms Mf are 
shown in Fig. [ 5 ] As mentioned above, there are two 
possible placements for the two SOC legs, in opposite 
sides of the loop. For each diagram, one has to also con¬ 
sider its hermitian-conjugate partner, since the interac¬ 
tion vertices and the SDW vertices are not adjoint oper¬ 
ators. Furthermore, the diagrammatic rules derived for 
this problem impose an overall minus sign to each one- 
loop diagram, and enforce the trace over the Pauli ma¬ 
trices to be taken in the direction opposite to the arrows. 
Finally, to compute these traces, it is useful to employ 
the following Pauli matrix identity: 


tr (ctVW) = 2 (S ij 6 kl - 6 ik 6 jl + 6 il S jk ) (31) 

A straightforward evaluation of these four diagrams 
gives then three different anisotropic terms for the mag¬ 
netic free energy: 
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FIG. 5. Illustration of distinct second-order one-loop diagrams in the case where the electron-electron interaction is given by 
V 13 and the magnetic order parameter is Mi. Note that the electron-electron vertex depends on the direction of momentum, 
i.e. the upper left and upper right diagrams are not identical. 


5F=- 


visy ^ ( tr (T 1 G' x T-G' T T'>G' T f 1 G X T-G r T''Gr) + tr ^T 1 G , r T y G , r T + G' x f 1 GrT y GrT + Gx) ) (M 2 2 - M 2 X - M 2 J 

kk' 

OJ n U) n / 

Vl3 Y ^ ( tr (T 1 G' X T-G' Y f 3 G’ r T 1 GxT-G Y f 3 Gr) + tr (t 1 G' T YG' y t + G' x r v G r r 3 G Y T + G x ) ) (Mf x - M\ y - M 2 Z ) 

kk' 

OJ n Ld n / 

visy ^ ( tr (^G'xT+G'y^G^Gxt+Gy^Gt) + tr {f 1 G' r f 3 G' y t~G' x f 3 Gyt~G x ) ) (Mf y - M? >x - M? iX ) 


kk' 

UJn UJ 


SF = u 13 A 2 ( ^ [Gx]i2[Gy]i2[Gt]2i) 2 (M 2 x - Ml y - M- 


U ■ 


(32) 


k,ct;„ 


Here, the primed Green’s functions are short-handed 
notations for G' = Gi (k',u4). The Green functions are 
given by 


lr i a m( k )^ a m( k )A 

[ AjiJ ^ iuj n - (e Aro (k) - ft) 


(33) 


where eA,m{ k) are the eigenenergies of the matrix }\a (k) 
and a l m (\i)A is the unitary transformation between the 
spinor basis and the band basis. In Eq. ( |32| ), the coeffi¬ 
cient of the term (M'f z — — Ml ) vanishes because 

of the antisymmetry of the GrT^Gr matrix in spinor 


space, whereas the coefficient of (M'l y — Mf — Mf z ) 
vanishes because [Gr] 12 (k) oc k x k y , causing the sum over 
momentum to vanish. 

Note that while the upper diagrams shown in Fig. [5] 
introduce the anisotropy between the in-plane and out- 
of-plane components of the magnetic order parameter, 
the lower diagrams contribute to the anisotropy between 
Mi iX and M-\ y . The reason for this is the character of 
the SOC in the effective model, which remains diagonal in 
the spin sector near the T-point due to xz- and yz -orbital 
characters of the electronic states. As a result, particle- 















hole excitations involving fermions from the T—point and 
from the X/Y -point do not allow for a spin-flip, which 
would be necessary to generate the anisotropy between 
M-| x and M-\ y . Only the inclusion of the particle-hole ex¬ 
citations between the two electron pockets, as described 
by the lower diagrams, generates the anisotropy between 
the x and y components of the magnetization. 

Repeating the same calculation for M 2 gives the same 
result, but with M\ >x —> M 2iV and M\^ y —> M 2 , x , as 
expected by symmetry. Therefore, we can recast this 
contribution to the free energy in the form of Eq. (30) 
via differences in the anisotropic coefficients a: 


„(«13) 


„(ri 3 ) 


-a — - 


a'” 13) - c4* 13) = 0 . 


2 ^i3A|f ( ^ [GA']i2[Gv]i2[Gr]2i(^84) 

(35) 


k,a; n 


The same procedure applied to the other interactions V 15 , 
i’n, and tqg reveals that only the first two give rise to 
anisotropic terms. The final result for the anisotropic 
coefficients is: 

a 2 — a.\ = 2t!i 5 A 2 ^ ^ [Gr]n[Gx]n[Gy]22^ 

k ,oj n 

— 2?q 3 A 2 ^ ^ [Gr ]21 [Gx] 12 [Gy] 12 ) ,(36) 

k,C O n 

013 — c^i = 2vn \ 2 ^ ^ [Gr*T y Gr]i 2 [Gx]n^ 

k ,u) n 

— 2rq 3 A 2 ^ ^[Gr]2i[Gx]i2[Gy]i2^) -(37) 

k ,uj n 


Mapping these interactions back to the more familiar 
multi-orbital Hubbard model, as done in Ref. 135] . 
reveals that U 13 = V 15 = vn = J, while the non¬ 
contributing term is Vig = U/2. Thus, our micro¬ 
scopic calculation reveals that the low-energy magnetic 
anisotropy arises from a combination of the SOC and of 
the Hund’s rule coupling. This anisotropy is present in 
the paramagnetic tetragonal phase, and does not require 
orbital or nematic order. 

It is now straightforward to determine which of the 
three terms in the free energy Eq. (301 dominates. For 
instance, if both a 2 — ay > 0 and 03 — aq > 0 , aq is the 
smallest of the three coefficients and the ordered com¬ 
ponents of the magnetic moments will be and/or 

M 2 y . By evaluating the expressions in Eqs. (36) and 
@ numerically, using the parameters that give the band 
dispersions and Fermi surface of Fig. [2] we can establish 
an effective “doping-temperature phase diagram” for the 
dominant anisotropy term as function of different values 
of the chemical potential /r and of the magnetic transi¬ 
tion temperature T. Because the phase boundaries of 
this phase diagram are given by the conditions a 2 = a 1 
or 03 = a 2l and because these coefficients are indepen¬ 
dent of U and have J as an overall pre-factor (see Eqs. 
(36) and (37)), the phase boundaries do not change by 


varying U and J. The phase diagram, shown in fig. [ 6 j 
reveals that for most of the parameter space considered 
here, the aq term is the smallest one, implying that the 
magnetic moments point parallel to their ordering vec¬ 
tors Qi below the magnetic transition. There is a small 
range of parameters in which the moments lie in-plane, 
but perpendicular to their ordering vectors (i.e. a 2 is 
the smallest). Such a parameter regime is likely not rele¬ 
vant for the iron pnictides, since it would require an “un¬ 
doped” composition (i.e. y, = 0 ) to display a rather small 
magnetic transition temperature. Most interestingly, we 
find a robust region in which the moments point out- 
of-plane (i.e. 03 is the smallest). This happens at any 
temperature, but always in the hole-doped side of the 
phase diagram (fj, < 0). In Fig. [7J we plot a zoom of the 
behavior of a 2 — aq and 03 — oq as function of the chem¬ 
ical potential for two fixed temperatures to illustrate the 
different regimes obtained. Note that, for most of the 
phase diagram, ( 0:2 — oq) <C ( 0:3 — oq), regardless of the 
value of J A 2 , which appears as an overall prefactor of all 
aq terms. This implies that the spin anisotropy behaves 
effectively as an easy-plane anisotropy. The fact that 03 
becomes the smallest coefficient in a narrow region of 
the phase diagram can be attributed to the fact that the 
term Z)k,a;„[ G! ri'r 2 / Gr]i 2 [Gx]ii in Eq. |37| changes sign 
from hole-doping to electron-doping. This behavior can 
be understood qualitatively by considering a hypotethical 
band structure in which all pockets are perfectly nested, 
Gr = (ico n + e) ® t° and Gx = Gy = (*w n — e) <g> r°. A 
straightforward calculation reveals that the three-Green’s 
function term above has different signs for y > 0 and 
/r < 0 , implying that it must vanish for a certain chemical 
potential value in the case of a realistic band structure. 

Our results reveal not only an important asymmetry 
between electron- and hole-doping, but also agree quali¬ 
tatively with experiments. In particular, neutron scatter¬ 
ing measurements in hole- and electron-doped BaFe 2 As 2 
f35U38| find generally in-plane moments parallel to Q, in 
the magnetically-ordered state, except at a narrow hole- 
doping range in which the moments reorient and point 
out-of-the-plane. 


IV. TETRAGONAL VS STRIPE MAGNETIC 
ORDER 


The previous section established the direction of the 
magnetic moments, but not the magnetic ground state. 
For instance, from the analysis of the second-order terms 
of the free energy, it is impossible to distinguish the cases 
in which either Mi or M 2 condense (i.e. M 2 = 0 or 
Mi = 0 ) from the case in which both condense simul¬ 
taneously (Mi = M 2 7 ^ 0). The former case gives the 
striped orthorhombic magnetic phases shown in Fig. [3j 
whereas the latter case gives rise to a double-Q (i.e. bi¬ 
axial) magnetic state that preserves the tetragonal sym¬ 
metry of the system. From the form of the anisotropic 
terms in the free energy (see Sec. Ill), there are three 
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FIG. 6 . Doping-temperature phase diagram displaying the 
smallest magnetic anisotropy coefficient a. When ai is the 
smallest, the moments point in-plane and parallel to the or¬ 
dering vectors; when Q 2 is the smallest, the moments point 
in-plane but perpendicular to the ordering vectors; finally, 
when Q 3 is the smallest, the moments point out-of-plane. The 
corresponding uniaxial configurations are shown in Fig. [3] 
Note that temperature here actually refers to the magnetic 
transition temperature, as our model approaches the onset of 
long-range magnetic order from the paramagnetic state. 


Act 



^ (meV) 

FIG. 7. Plots showing the dependence of 03 — cm (red 
curve) and 02 — 01 (green curve) for two distinct temperatures 
(see the dashed lines in Fig. pj as function of the chemical 
potential. 


different types of tetragonal magnetic ground states, as 
shown in Fig. [ 8 j Two of them correspond to the so-called 
spin-vortex crystal phase (SVC), a non-collinear state in 
which = M 2 tV 7 ^ 0 or M\^ y = M 2 , x 7 ^ 0 , whereas 

the third one corresponds to the so-called charge-spin 
density-wave phase (CSDW) [ST], a non-uniform state in 
which Mi j2 = M 2i2 7 ^ 0 . 

To determine whether the ground state corresponds 
to an orthorhombic uniaxial SDW (i.e. either Mi 7 ^ 0 


SVC CSDW 

OL 1 OL2 (A3 


\ / \ 

/ \ 

/ \ 


• 0 • 

/ \ / 

\ / 

\ / 

\ * 

0 « 0 

\ / \ 

/ \ 

/ \ 

✓ ® 

. 0 • 


/\/\/\/ \ • ® ® 

FIG. 8 . Sketch of the three possible biaxial tetragonal mag¬ 
netic phases. Here SVC is the spin vortex crystal phase 
(with non-collinear magnetic moments) and CSDW is the 
charge-spin density wave phase (with non-uniform magnetic 
moments). 


or M 2 / 0) or to a tetragonal biaxial SDW (i.e. Mi = 
M 2 7 ^ 0), we need to go to higher order in the free energy 
expansion. Symmetry requires the free energy to have the 
form (in the absence of SOC): 


F (4) (Mi, M 2 ) = | (Ml + M“) 2 - | (Mi - M|) ^38) 
+ 2io(Mi • M 2 ) 2 . (39) 


Minimizing this expression shows that the tetragonal 
biaxial state is realized when g < 0 or g < —w, whereas 
the orthorhombic uniaxial state takes place when g > 0 
and g > —w. 

The same model was derived by different itinerant ap¬ 
proaches for the magnetic instabilities of the iron pnic- 
tides, revealing different parameter regimes in which 
the uniaxial or the biaxial states are the ground states 
[El EU [22l HU [52H54]. In this regard, the novelty of our 
approach relies on the relationship between these ground 
states and the magnetic anisotropies, and also on the em¬ 
ployment of a low-energy model that respects all symme¬ 
tries of the FeAs plane, including the As puckering that 
enhances the size of the Fe unit cell. Because u and g 
are non-zero even for vanishing SOC and interactions, 
we compute only the contributions arising from the non¬ 
interacting part of the Hamiltonian. This is achieved 
either by standard diagrammatics or by explicitly inte¬ 
grating out the electronic degrees of freedom. We find: 
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FIG. 9. Illustration of the three distinct fourth order diagrams that contribute to F^ 4 K Note that the right diagram has a 
symmetry factor of two, while the symmetry factor of the middle and left diagrams is one. One might expect there to be a 
fourth diagram with alternating Mi and M 2 , however, the vertex coupling the SDW order parameter and the fermions forbids 
it (see Appendix [P] for more details). 


E (4) = M 4 J2 tr (^Grf'Gx^Gr^Gx) 

k.c j n 

+Mj Y tr (r+Grr-Gyr+Grr-Gy) 

k ,uj n 

+2MJM2 J2 tr (' r+ GrT 1 GxT 1 G r r _ Gy) 

+M2 Y, ([Gr ]| 2 [^V]ll) 

k ,uj n 

+ 2 M?M^ Y ([Gr] 12[Gx] 11 [Gy] n ) , ( 40 ) 

k ,u) n 


with the corresponding diagrams shown in Fig. [9] 
Rewriting the free energy in the form (381, we can readily 
obtain u, g , and w : 


u - (I^r] 12 [Gy] 11 [Gx] 11 + [Gr}u[Gx]li) (41) 

k ,OJ n 

9 = Y ([^r]i 2 [G f y]ii[Gx]ii - [Gr]?i[Gx]?i) (42) 

k.cD^ 

w = 0. (43) 


One might expect the third diagram in Fig. [9]to result in 
a non-vanishing w, however, contraction of the Pauli ma¬ 
trices (see Eq. 31) reveals that the Mi • M 2 term cancels. 
More generally, w = 0 is a robust property of our model, a 
consequence of momentum conservation and the absence 
of a Fermi pocket at (n,n) (see Appendix [dJ. Including 
a hole-pocket at ( 7 r, 7r) or interactions will however lead 
to a non-zero contribution to w |. r >31 . 

We can now compute numerically the value of g for 
the same “doping-temperature phase-diagram” studied 
in the previous section, see Fig. [6] The combined re¬ 
sult, shown in Fig. [lO] accentuates the asymmetry be¬ 
tween hole- and electron-doping discussed previously. In 
particular, while hole-doping tends to favor a tetrago¬ 
nal biaxial SDW state, electron-doping tends to favor an 


orthorhombic uniaxial SDW state. A cut with the be¬ 
havior of g as function of g for a fixed temperature is 
also shown in Fig. 11 To gain more insight into the 
behavior of g , we consider once again the hypothetical 
case of perfectly-nested bands, Gy = (*w„ + e) g> r° and 
Gx = Gy = ( iujn — s) <8S T°. In this case, from the equa¬ 
tions above, g < 0. Building on the results of Ref. US], 
we expect that the sign of g will change once the two 
hole pockets become rather different in size, such that 
one of them becomes poorly nested with the electron 
pockets. Our calculations indicate that, for the general 
tight-binding model studied here, this is favored by hole 
doping rather than electron doping. 

It is important to emphasize that these results should 
be understood as general trends as function of the chem- 



stripes /i (meV) 

FIG. 10. Doping-temperature phase diagram of the dif¬ 
ferent types of magnetic ground state (stripes or tetragonal) 
and their corresponding spin orientation (in-plane or out-of- 
plane). The color-code corresponds to the magnetic configu¬ 
rations shown in Figs. [3] and [8] Note that temperature here 
actually refers to the magnetic transition temperature, as our 
model approaches the onset of long-range magnetic order from 
the paramagnetic state. 
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FIG. 11. The quartic coefficient g as a function of the chem¬ 
ical potential g for a constant temperature (see the dashed 
line in Fig. 101. When g is negative the system chooses a 
tetragonal biaxial magnetic phase whereas for positive g, the 
system selects an orthorhombic uniaxial stripe state. 


ical potential, rather than a full determination of the 
ground state for each specific value of g. This is because, 
in contrast to the previous section, in which the lowest or¬ 
der contribution to the spin anisotropy arises solely from 
the SOC and the interactions, there are other potential 
contributions to g beyond the scope of the current work. 
Among these contributions, we highlight the sizable 
magneto-elastic coupling, which should extend the stripe 
phase to wider doping ranges { 55 ] [ 56 ] , and interaction cor¬ 
rections, which can also favor the uniaxial over the biaxial 
state mm- With this word of caution, we note that the 
tendency observed here that hole-doped compounds are 
more favorable to a tetragonal magnetic phase as com¬ 
pared to their electron-doped counterparts is in qualita¬ 
tive agreement with experiments, which observe a small 
region of tetragonal magnetism near the optimally-hole 
doped pnictides Bai_ x Na a; Fe 2 As 2 Ml Sri_ 3 ; Na a; Fe 2 As 2 

[35] , Bai_ x K x Fe 2 As 2 {55} [51], and Ba(Fei_ a ,Mn ;c ) 2 As 2 

EDI- 


Interestingly, in the first three compounds, neutron 
scattering has shown that the onset of a tetragonal mag¬ 
netic state takes place in a region in which the mag¬ 
netic moments reorient from in-plane to out-of-plane. 
Within our analysis, this can be attributed to the robust 
region in parameter space in which 03 is the smallest 
anisotropic coefficient (see Fig. [lO]). More importantly, 
this anisotropic coefficient removes the degeneracy be¬ 
tween the two types of tetragonal SDW phase - the SVC 
and the CSDW states (see Fig. [ 8 ]) by favoring the lat¬ 
ter. This is expected to happen even if w ^ 0 , since the 
latter is a quartic coefficient, whereas 0:3 is a quadratic 
coefficient. Therefore, at least near the onset of the mag¬ 
netic transition, it is the SOC and the Hund’s rule cou¬ 
pling that select the CSDW phase. Recently, Mossbauer 
[351 and p,SR [5T. experiments in Sri_ x Naa;Fe 2 As 2 and 
Bai_ a ,K ;c Fe 2 As 2 , respectively, have reported direct ev¬ 
idence that indeed the CSDW state is realized in the 
regime where the spin is reoriented and the magnetic 
long-range order preserves tetragonal symmetry. 


V. DISCUSSION AND CONCLUSIONS 


In summary, we have shown that within a low-energy 
model that respects the symmetries of the FeAs plane, 
magnetic anisotropy arises naturally from the combina¬ 
tion of the spin-orbit coupling and the Hund’s rule cou¬ 
pling. The magnetic anisotropy consists of three terms 
(see Sec. Ill I. Although it cannot be mapped generally on 
an easy-axis or an easy-plane term, it effectively behaves 
as an easy-plane term for a large part of the parameter 
region studied here, since (02 — aq) -C ( 0:3 — 01 ). We 
found that, for most of the temperature-doping phase 
diagram, the spin anisotropy is such that the magnetic 
moments point in-plane and parallel to the direction 
of the ordering vector. For a small doping range in 
the hole-doped side, across all temperatures studied, 
the magnetic moments tend to reorient and point out- 
of-plane. These features are consistent with those ob¬ 
served experimentally, including the spin reorientation 
observed in the hole-doped pnictides Bai_ a; Na a: Fe 2 As 2 
[35] . Sri_ x Naa;Fe 2 As 2 [35], Bai_ a; K x Fe 2 As 2 (37]. We also 
found a general tendency of tetragonal double-Q mag¬ 
netic order for the hole-doped side of the phase diagram, 
whereas the orthorhombic single-Q stripe magnetic or¬ 
der is favored in the electron-doped side. Although this 
is in general agreement with the experimental obser¬ 
vations in Ba 1 _ a ,Na a; Fe 2 As 2 [57], Sr 1 _ a; Na a; Fe 2 As 2 [35] , 
Bai_ a ,K x Fe 2 As 2 [55], and Ba(Fei_ a; Mn 2 ,) 2 As 2 [50] , our 
results seem to overestimate the size of the region in 
which the tetragonal magnetic state is stable. A pos¬ 
sible reason for this discrepancy is that our model does 
not account for other factors that usually favor the stripe 
over the tetragonal magnetic state, such as the magneto¬ 
elastic coupling [55] [55] and the residual interactions not 
directly responsible for the SDW instability mm- Yet, 
our results provide a clear connection between the spin 
reorientation and the type of tetragonal magnetic state 
observed in the hole-doped iron pnictides - namely, the 
charge-spin density-wave state with a non-uniform mag¬ 
netization OH El]. 


An important consequence of our results is that the 
spin anisotropy is not necessarily tied to the orbital order 
that is triggered across the nematic/structural transition 
[ 30 ] . Although it is plausible that such an orbital order 
affects the spin anisotropy, the latter exists already in the 
tetragonal paramagnetic state as a result of the symme¬ 
try properties of the FeAs plane. In this regard, it would 
be interesting to investigate how the spin anisotropies in 
the tetragonal-paramagnetic phase studied here are con¬ 
nected to the spin anisotropies in the low-temperature 
phase, after both magnetic and nematic orders are well 
established. Finally, our results open the important ques¬ 
tion of how this particular form of magnetic anisotropy 
impacts the normal state properties of the iron pnictides. 
In particular, the onset temperatures and the characters 
of the coupled nematic-magnetic transitions are expected 
to be strongly affected by any form of spin anisotropy 

m ih iMj- 
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Appendix A: Expansion of tight-binding 
Hamiltonian for small k 

In this appendix we derive explicitly the non¬ 
interacting Hamiltonian Hq introduced in sec. [TlJ As 
explained in that section, we need to project and ex¬ 
pand the 5-orbital tight-binding dispersion £ yi/ (k), where 
fx, v = 1,..., 5 are the orbital indices corresponding to xz , 
yz, x 2 — y 2 , xy and 3z 2 — r 2 , respectively. The dispersions 
are given by [46] 


£n = txz/yz + 2f 41 cos k x + 2f 41 cos ky + 4f 44 cos k x cos k y + 2(cos 2 k x - cos 2 k y ) 
+ l cxy cos 2k x cos k y + 4f 44 y cos k x cos 2 k y + 4f 44 cos 2 k x cos 2 k y , 

£22 = Cxz/yz + 2f 22 cos k x + 2f 22 cos k y + 4f 22 cos k x cos ky - 2t 2 xx / yy (cos 2 k x - cos 2ky) 
+ 4 t 22 xy cos 2k x cos k y + 4 t xyy cos k x cos 2 k y + 4f 22 yy cos 2 k x cos 2 k y , 

£33 = £x 2 -y 2 + 2 t 33 / y (cos k x + cos k y ) + U 33 y cos k x cos k y 2t 33 x / yy (cos 2 k x + cos 2 k y ) , 
£44 = Cxy + 2f] l f/ y (cos k x + cos ky) + t xy cos k x cos k y + 2 t xx , yy (cos 2 k x + cos 2 k y ) 

+ 4 t xxy / xyy (cos 2 k x cos k y + cos k x cos 2 k y ) + 4t xxyy cos 2 k x cos 2k y , 


£55 

— Ez 2 4 

- 2t l% 

(cos k x 

+ cos , 

ky) + 2t 

ll/yy (cos 2k x cos2k y ) 




+ Kly/xyy (COS 2k x COS k y ■ 

■f cos k x 

cos 2k y ) + 4 :t 5 x 5 xyy cos 2 k x cos 2 k y , 


£12 

= <' 

sin k x i 

sin k y + 

4 tlly/xyy ( sin 

2 k x sin k y + sin k x sin 2 k y ) 

1 -L 4 t 12 

1 xxyy 

sin 2 k x sin 2 k y 

£13 

= z2tf 

sin k y 

+ 

cos k x 

sin k y - 

-iKly/xyy (COS k x sin 2k y 

— cos 2 k x 

sin k y ) , 

£14 

= i^x 

sin k x 

~ iKt 

sin k x 

cos k y + sin 2 k x cos k y , 



£15 

= i2tf 

sin k y 

- 

cos k x 

sin k y - 

- iAtllyy cos 2 k x sin 2 k y , 



£23 

= *2 t 23 

sin k x 

+ 

sin k x 

cos k y - 

- iAt lly/xyy ( sin 2k * COS k V 

— sin k x cos 2k y ) , 


£24 = ~®2f 24 sin k y + i4f 24 cos k x sin k y — *4f 24 y cos k x sin 2 k y , 

£25 = — *2t 25 sin k x + i4:t xy sin k x cos k y + i4:t x 5 xyy sin 2 k x cos 2 k y , 
e 34 = 4 tl xy / xyy (sin k x sin 2 k y - sin 2k x sin k y ) , 

£ 35 = 2 t 3 x/y (cos k x - cos k y ) + 4 t 3 xxy/xyy (cos 2k x cos k y - cos k x cos 2 k y ) , 
£45 = 4 t xy sin k x sin k y + 4 t x 5 xyy sin 2 k x sin 2 k y . 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(AS) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 
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Here e* are the onsite energies associated with each orbital and are hopping parameters from orbital fi on site i 
to orbital v on site j. The above expressions are accompanied by constraints on the coefficients tY' due to tetragonal 
symmetry: 


** = * 
HI 


t 11 = t z 

V 

j.11 __ y-22 

^xyy xxy 
£13 — +23 

L xxy/xyy h xxy/xyy 
+15 _ f2b 
L y ~ L x 


+22 


4.11 _ fZZ 

^xx^ L xyy 

^xy ^xy 


+24 


t = t 

u xxy xyy 


t 11 = t Z 

xy xy 
4-11 _ 4.22 

''xx'L^ L xxyy 

t X = ty 

+15 _ +25 
L xy L xy 


+22 


4.11 _ 4.22 ,xx/yy 

xx/yy 

4-13 _ 423 

Ly L X 

4.14 _ 4.24 

+15 _ +25 

L xxyy L xxyy * 


(A16) 


We are now in a position to expand the elements of e M „(k) around the T, X and Y points. At the T point the 
orbitals xz and yz dominate, corresponding to the elements £n, £ 12 , £21 and £ 22 - Similarly, at the X ( Y ) point the 
dominant orbitals are yz (xz) and xy. To obtain these parts we expand £22 (£n), £24 (£ 14 ) and £44 around (k x + n, k y ) 
-tt)): 


(( k x ,k y 


hr = 


h x = 


hy = 


C 1 + C 2 (k 2 x + k 2 y )+C 3 (k 2 x -kl) 

C/^kx ky 

C 5 +C 6 (k 2 x + k 2 y )+C 7 (k 2 x -k 2 y ) 

iv x(k) 

_ fC 5 + C 6 (k 2 x + kl)-C 7 (k 2 -k 2 y ) 


ivy( k) 


C A k x k v 

C 1+ C 2 (kl + k 2 y )-C 3 (k 2 ^k 2 y ) 
-iv x ( k) 

C 12 (k 2 x + k 2 y )+C 13 (k 2 x -k 2 y ) 

—iuv(k) 

C 12 (k 2 x + k 2 y ) - C 13 (hi - k 2 ) 


C 


11 


c 


where 


«x(k) = C 8 k y + C 9 k y ( k 2 y + 3 k 2 x ) - C w k y (k 2 x - k 2 y ) 


uv(k) = ~C 8 k x - C 9 k x ( k 2 x + 3 k^) - C w k x ( k 2 x - k%) 

As a function of the tight-binding parameters, the constants C 1 ,..., C 33 are 

Ci = er = txz/yz + 2 (tl 1 + i/) + 4 ( t]/ xy + t X xyy + t]) y + t l / xyy ) 
C 2 = 2 


-^ = -\ (*“ + t l 1 )- 5 (/xxy + tlly) - 2 tH - 8t - 


11 

xxyy 


C 3 = b = \ l 1 - 4 1 ) + 3 (Cj/ - *xxy) - 4i : 


11 

xx/yy 


_ a „ _ _ a f 4.12 1 /I+12 1 4+12 \ 

'-'4 ^ l ^xy ' ^xxyy ^ xxy/xyy I 


C 5 = Cl = e-xz/yz + 2 (4 1 - iy 1 ) - 4 (4xy ~ 4yy + 4i xy “ 4i 

C 8 =2 


11 ) 
xxyy) 


(*“ - *“) + 5 (t“y - 4yy) + 2 C “ 8t 


11 

xxyy 


1 


C 7 = ai = - (tl 1 + t 1 /) - 3 (tHy + illy) + 4t . 

C 8 = 2v = 2 (tj , 4 + 2f 44 - 4t 44 y ) 

_j_+ 24 _ ^424 17 

12 v 6 

Cio = 2 P2 = — zt 


11 

xx/yy 


424 


C g = 2pi = ~—t y - —t x l + -t z xyy 


y.24 1 q+24 
4 l V ' OL xyy 


C \1 — e 3 — e xy + 4 (-t 4 ^ + tii/yy + iiiyy'j 

1 t .AA ~.AA - 1 44 \ 

V 


C12 = 2 


— 0 (+44 _ 2/44 _ 444 

2 m 3 \ xy xx /w xxyy ; 


444 


1^13 ^3 ix/y ^ixxy/xyy ' 


(A17) 

(A18) 

(A19) 

(A20) 

(A21) 

(A22) 

(A23) 

(A24) 

(A25) 

(A26) 

(A27) 

(A28) 

(A29) 

(A30) 

(A31) 

(A32) 

(A33) 

(A34) 


The overall minus sign in the coefficient C 4 arises due to the minus sign in the definition of the spinor in Eq. Q . The 
coefficients can be obtained either by using the relations above with the coefficients determined from tight-binding 
fits to DFT calculations, or by directly fitting the coefficients to DFT calculations. 


Appendix B: Spin-orbit coupling in orbital basis 


of the main text. Denote the eigenstates of L z by \m) 


Here we express the standard spin-orbit coupling term 
AS • L in the orbital basis, which leads to Eqs. (19)-(21) 
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where m = The spin-orbit Hamiltonian is 

then 

Hsoc = ^{ma |AS • L| nP)<f ma d n p . (Bl) 

mn 

a(3 

where d' ma creates an electron with spin a and angular 
momentum projection to. Using the fact that S ■ L = 
L Z S Z + \{L+S- + L~S+) and 

S z \na) = ±-|na) , S' ± |na) = \n, a ± 1) ,(B2) 

the Hamiltonian becomes 

Hsoc = [( m \ L z\ n ) d ln,t d ™t - {™\ L z \ n )d ] mi d ni 

mn 

+ (m\L + \n)d' mi d , nt + {m\L~\n)d' m ^d n i 
= ^Y, A ™ndl a dn0, (B3) 

mn 

a/3 

with the matrix elements: 

A mn = ~ A mn = ( ™\ L A n ) = n5 mu (B4) 

A rnn = \/ {L — n)(L + n + l)<5 m ,n+l (B5) 

A mn = 7(-^ + n )(L — Tl+ l)S m ,n-l ■ (B6) 

We can transform the Hamiltonian in Eq. |B3| to the basis 
spanned by the cubic harmonics (i.e. the orbital basis) 
using 


^ma — ^ ^ Um^C^Q^ 
n 

(B7) 

S = V U* c f 

ma / j ^ mfi^fia 5 

(B8) 

where U m ^ = (m\y) and \y) is the basis states in the 
space of cubic harmonics. The transformed Hamiltonian 
is 

Hsoc = ^ ^2 A^c( a c v p , 

fils 

a/3 

(B9) 

with A a/3 = WA^U. 

To proceed we specialize to the case L = 2, 
in the well-known d-orbitals. The basis states 
(xz\, {yz |, (xy |, (x 2 -y 2 \,(z 2 \ with 

resulting 
are (y\ = 

( TO | = i=(-(l| + <-l|) 

(B10) 

( yz \ = ^(-(l\-(-l\) 

(BH) 

( X y\ = ^(-(-2\ + (2\) 

(B12) 

(x 2 -y 2 \ = ^=((-2\ + (2\) 

(B13) 

7 2 l = (0|, 

(B14) 


and W can be read off as the coefficients of these equa¬ 
tions. Thus, we find U to be 


U=^ 

72 


/ 0 0-il 0 \ 
-1 i 0 0 O' 
0 0 0 0 ^2 
1*000 
V o o * 1 o / 


(B15) 


yielding the three independent H-matrices: 


Itt = 


A* = 


Itl = 


/0 —* 
* 0 

0 

0 

0 

0 

0 




0 0 

0 

*2 

0 

5 


(B16) 

0 0 

—*2 

0 

0 




\0 0 

0 

0 

0 ) 




(0 

0 * 

1 


-73\ 



0 

0 1 

— 


-*73 



—i 

-1 0 

0 


0 

1 

(B17) 

-1 

* 0 

0 


0 



\V3 * 

73 0 

0 


0 ) 



( 0 

0 

* 


1 7s\ 


0 

0 

-1 


-i —i 

73 


—i 

1 

0 

0 0 

. (BIB) 

1 

i 

0 

0 0 


\-73 

*73 

0 

0 0 ) 








(B19) 


Considering only the \xz), \yz) and \xy) orbitals, cor¬ 
responding to the upper left 3x3 blocks in the above 

-— a/3 

matrices results in the restricted matrix A '^, v ,: 

/ 0 -i 0 0 0 * \ 
i 0000-1 
0 0 0 -i 1 0 

0 0 * 0*0 
0 0 1 -* 0 0 
V-* -i o o o o / 


A! = 


(B20) 


Decomposing it into spin and orbital parts gives: 
A' = 


+ 


0 

0 

M 

1 

(o 

0 

0 

0 

0 

0 

<S> <T X + | 

0 

0 

—* 

—* 

0 


f 

v° 

* 

0 

0 

—* 

o\ 

i 




* 

0 

0 

\ ®<T Z . 




0 

0 

0 

1 





(B21) 


Finally, applying this expression to the Hamiltonian 


(|B9j) gives 

sU 


M V U /X- 


c v ‘ — * ^ c! cz cl a a pC X y ! i 3 + h.c. 


.A 


+i~ct„ „al a c vz . B + h.c. 


' 2 l,a u 

X 

c 

2 


+7 C yz,a a a0 C xz,0 + ^.C. ,(B22) 


which leads to Eqs. (19) (21). 
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Appendix C: Evaluation of Two-loop Diagrams 


Diagrams of the two-loop type, as shown in Fig. Bh 
can be split into irreducible diagrams. As a result, all 
two-loop diagrams for a given interaction can be ob¬ 
tained by squaring the sum of irreducible diagrams, as 
illustrated in Fig. [12] 

To illustrate the cancelation of the two diagrams in 
the sum, we choose the interaction iq 7 and the order 
parameter Mi. In this case, the sum in the brackets in 


Fig. 12 is, for the SOC leg related to the T spinor: 


A 


- ^2 <j l ]tv[T + G tt v Gtt 1 Gx] 

k,n 

M i,ite[< jZ(jl ]to[ T ~ G xT 1 G T T v G T ] 

k,n 

=A^M M ([G r r ! 'Gr]2i[G 

n 

+ [Gyt v Gr]i2[Gx]n) 5 


*Xlll 


k.n 


(Cl) 


This term vanishes since GrT y Gr is an antisymmetric 
matrix. For the contribution from the diagrams with the 
SOC leg related to the X/Y spinors, we find: 

- *A ^2 (Mi iX tr (f 1 Gvr _ Gj f r 1 Gr) 

k ,n 

+M ltV tr (f'GYT+Gxf'Gr) ) 

= —*A^ (■^i,x[Gv]i2[Gx]ii[G f r]n 

k,n 

+Mi i j / [Gv]n[Gx] 2 i[Gr]n^ , (C2) 


which is also zero as the off-diagonal elements of the 
Green functions Gx and Gy are odd functions of k. 
Similar arguments apply in the case when the electron- 
electron interaction is given by either V 13 , V 15 or ^19 and 
when the magnetic order parameter is M 2 . Thus, all con¬ 
tributions from the two-loop diagrams vanish, and we are 
left with only the one-loop diagrams shown in Fig. [4[>. 


Appendix D: Diagrams contributing to w 


In this appendix we explain in more details the state¬ 
ment made in the main text concerning the vanishing of 
the quartic coefficient w in Eq. (38 1 . Let us consider a 



FIG. 12. Illustration of the decomposition of a two-loop 
diagram into irreducible diagrams. 


Mi • a 


M2 • 0 


0 . 0 ) 


G(k + Qi) > t 


G(k + Q 2 ) 



FIG. 13. Illustration of a generic diagram contributing to 
the coefficient w. Here Qi = (n, 0) and Q 2 = (0,7r). On the 
right diagram we have imposed momentum conservation at 
each vertex, resulting in the appearance of G(k + Qi + Q 2 ). 


generic diagram which would contribute to the coefficient 
w , as shown in Fig. [13] 

Note that we do not specify any vertices, i.e. the cou¬ 
pling between the electrons and the SDW order param¬ 
eters do not necessarily arise from Eq. [28] Due to the 
Pauli matrix contraction, the diagram must have alter¬ 
nating Mi and M 2 legs in order for it to contribute to 
w. Indeed, performing the trace over spin indices gives: 

M{M J 2 MiM l 2 tr (oVcrV) . (Dl) 

Using Eq. [31] we find: 

2 (M x -M 2 ) 2 -MjM?, (D2) 

thus resulting in a (Mi • M 2 ) 2 term. This contrasts to 
the third diagram in Fig. J9] which gives no contribu¬ 
tion of the form (Mi • M 2 ) after tracing over the Pauli 
matrices. 

Let us now consider the internal lines of the diagram. 
Since Mi carries momentum (n, 0) and M 2 carries mo¬ 
mentum ( 0 , 7 r), the only way for momentum to be con¬ 
served is if one of the lines corresponds to a propaga¬ 
tor with momentum Qi + Q 2 = ( 7 r, 7 r). However, in 
the absence of a Fermi pocket at M = (n, n) (of the 
unfolded Brillouin zone), this will be an off-shell contri¬ 
bution. Thus, contributions to w must arise from the 
electronic states near M = ( 7 r, 7 r). 















